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SUMMARY 

This paper presents a closed-form expression for the weight function for 
a strip with a single edge crack. The expression Is valid for relative crack 
lengths fr«n zero to unity. It Is based on the assumption that the shape of 
an opened edge crack can be approximated by a conic section. The results agree 
well with published values for weight functions, stress Intensity factors, and 
crack mouth opening displacements. 


INTRODUCTION 

Bueckner's weight function concept (refs. 1 and 2 ) has proven to be a use- 
ful tool for elastic crack analysis. The weight function Is a function of spe- 
cimen (or structure) and crac>. onfiguratlons only and Is Independent of the 
load system. If an expression for the weight function Is available, the stress 
intensity factor for any load system may be computed from the stress distribu- 
tion in the uncracked body by simple Integration. In presenting this concept 
Bueckner used the example of a strip with a single edge crack, which Is a com- 
mon test specimen configuration and which also has application to real-world 
structural problwns. He presented (ref. 2) a simple expression for the weight 
function for this configuration which Is useful for relative crack lengths from 
zero to one-half. 

Rice (ref. 3) showed that weight functions can be calculated from crack 
surface displacements and stress intensity factors. Thus, If the stress Inten- 
sity factor and the complete crack displacement solution are available for a 
given cracked-body configuration under one loading condition (say, simple ten- 
sion), the weight function can be confuted. This configuration and loading Is 
usually called the reference solution. 

At about the same time, this author (ref. 4) proposed that the opening 
shape of an edge crack could be approximated by a conic section. In the pres- 
ent paper, the conic section approximation Is used to develop a closed-form 
expression for the weight function which is valid for all values of relative 
crack length, I.e., from zero to unity. 


ANALYSIS 

The generalized equation for the conic section (ref. 4) and the notation 
used are shown In figure 1. The conic section coefficients which are shown 
(m ■ -1, 0, 10, ») correspond respectively to an ellipse, a parabola, a hy- 
perbola, and a pair of straight lines. Using relations given In reference 4, 
the opening displacement (n)'may be written as 
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where E' is the effective elastic nradulus, Kj is the opening-mode stress in- 
tensity factor, and m is the conic section coefficient. Rice (ref. 3) showed 
that the weight function (h) may be written as 
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In comparing weight function formulations, note that Bueckner's form M(x) and 
Rice's form h{x) are related by 

M(x) IT h(x) 

Now a general expression for the weight function may be determined by differ- 
entiating equation (1). In doing so, note that the coordinate system in fig- 
ure 1 is located at the crack tip. If the crack extends an amount da, the 
coordinate system moves a like amount (i.e., dx = da). After performing the 
required differentiation, Bueckner's reduced form of the weight function (N(x)) 
may be written as 

. X fx /X dm m\ . ml . m /,x 

where x » a/W is the relative crack length and Y = Kj/a V? is the dimen- 
sionless stress intensity factor coefficient. Thus the weight function can be 
computed if we know how the conic section and stress intensity factor coeffi- 
cients vary with relative crack length for one form of loading, say, simple 
tension. 

Values of the conic section coefficient are compuated as follows. Equa- 
tion (5) of reference 4 gives the relationship between stress intensity factor, 
crack mouth displacement, and conic section coefficient as 

E'ng/aa - 2Y[(2 + m)/irl^^^ 



Thus if stress intensity and mouth displacement are known, the conic section 
coefficient can readily be determined. Since the tip radius of an opened crack 
is proportional to the square of the stress intensity factor, this approach is 
equivalent to fitting a conic section to the opened crack shape at the mouth 
and at the tip. 

In order to determine values of the conic section coefficient, values for 
the stress intensity and crack mouth displacement were taken from the litera- 
ture. Stallybrass (ref. 5) gives the stress intensity factor and crack opening 
displacement for an edge crack in a half-plane as 
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Y . 0.7930^^ - 1.988 

E’ng/aa « 0.7930 x - 4) - 2.908 

from which the conic section coefficient at x ■ 0 Is 

m - (8 - ir^)/(ir^ - 4) - -0.3185 

Keer and Freedman (ref. 6) give stress Intensity and opening displacement val- 
ues for various edge cracks In finite plates. However, their results at 
X ■ 0.85 appear to be Inaccurate and were Ignored (their value for the stress 
Intensity coefficient Is Y » 42.9, whereas Benthem and Kolter (ref. 7) give 

Y » 33.0). Benthem and Kolter (ref. 7) treat the edge-dam In a half-plane, 
which corresponds to the major limit of the single edge crack. From their 
work, the limit values (x « 1) of the conic section coefficient and Its slope 
are 

11mx_^lCm(l - x)l = 3.104 
11mj^^jd[m(l - x)]/dx - 8.034 

Values of the conic section coefficient were computed from the results of ref- 
erences ? to 7 using equation (5) of reference 4, which was given earlier. 

A polynomial 1n terms of relative crack length was derived from these 
coefficients as follows. First, the conic section coefficients were put into 
finite form and a polynomial fit by the method of least squares. Then the 
polynomial coefficients were adjusted by trial and error to provide the correct 
magnitude and slope at each extreme. The resulting polynomial Is 

m(l - X) . -0.3185 + 0.3185 x + 1.3765 x^ 6.2475 X^ -8.82 x^ + 4.3 x^ (4) 

The polynomial Is everywhere finite, has the correct magnitudes and slopes at 
the extremities, Is within 1 percent of the reference values, and Is shown in 
figure 2. 

Finally a wide-range expression for the stress Intensity coefficient for 
tension Is needed, and Benthem and Kolter's (ref. 6) should suffice. Since 
only tabular values are given In reference 6, the equation was rederived as 

Y - ^7(1 + 2x) (1 - 21.006 + 87.065 x + 

2 1 / 2 ] 

-(395.38 + 3618.5 x + 7584.88 x^) J (5) 

The term in brackets is the coefficient tabulated In reference 7 and agrees 
with the tabulations within 0.001. Now having wide-range expressions for the 
conic section and stress Intensity coefficients, the calculation of weight 
functions Is straightforward. 
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Once the weight function has been determined, the stress Intensity factor 
for any loading can be determined by simple Integration, as 

Kj ■ p(x/a) N(x/a) x”^^^ dx (6) 

where p(x/a) Is the stress distribution over the crack location In the un- 
cracked body. 


RESULTS AND DISCUSSION 

Weight function values computed from equations (3) to (5) and those given 
by Bueckner (refs. 1 and 2) are presented In table I. The conic-section weight 
functions are within 4 percent of Bueckner's for relative crack lengths of one- 
half and less. As a further check, stress Intensity coefficients for tension 
and pure bending were computed by Integrating the conic-section weight func- 
tions. These results are presented In table II, and are within 1.8 percent 
(tension) and 2.7 percent (bending) of those given by Benthem and Kolter 
(ref. 7) over the entire range of relative crack lengths from zero to unity. 
Crack mouth displacement coefficients are presented In table III and will be 
discussed later. 


An alternate approach to the approximation of weight functions was pre- 
sented by Petroski and Achenbach (ref. 8). They prescribe the crack opening 
displacement In the form 




1/2 ^ 1/2 ^ 


G(x) 




^ . 


(7) 


In which the term 6(x) Is determined by requiring that the stress Intensity 
factor, when determined by Integration as In equation (6), be Identical to the 
value used In equation (7). In other words, equations (6) and (7) must be 
self-similar. The approximate weight function Is then determined using equa- 
tion (2). This approach Is advantageous If crack mouth displacements are not 
available for the reference problem. However, if such displacements are avail- 
able, th 2 conic section approximation yields more accurate displacement coef- 
ficients, as might be expected. This can be seen in table III. The coeffi- 
cients computed from equation (1) are within 0.7 percent of Keer and Freedman's 
(ref. 6) up to x = 0.65 and within 2.1 percent at x » 0.80. As mentioned 
earlier, the results In (ref. 6) for x » 0.85 are probably inaccurate. For 
comparison, the Interpolation equations of Srawley and Gross (ref. <J) and Tada 
(ref. 10) give crack mouth displacement coefficients of 150.1 and 149.8, re- 
spectively, at X s 0.85. 


Even though these results may be satisfactory, some further attention 
should be paid to the matter of crack opening shape. Consider Stallybrass's 
solution (ref. 5) for the edge crack In a half-plane, which was mentioned 
earlier. Note that the conic section coefficient Is an exact value. Indepen- 
dent of his numerical computations. This might lead one to expect that the 
crack opening shape Is an exact conic section. However, this Is not the case, 
as can be seen In figure 3. Here the displacements from the conic section 


a 


model are compared with Bueckner's (ref. 1) and the earlier results of Wlggles* 
worth (ref. 11). Although they seem small in figure 3(a), the differences are 
significant, as can be seen in figure 3(b).* The conic-section displacements 
agree with Wigqlesworth's at the crack tip and mouth, but are larger In be- 
tween. They are smaller than Bueckner's at the tip and mouth but, again, are 
larger In between. Similar small but consistent differences can also be found 
in the weight functions themselves (see table I). The conic-section weight 
functions are almost everywhere larger than Bueckner's for x <,0.40 and 
smaller for x £ 0.45. These small unidirectional errors do n'ot average out 
when the weight function is integrated to determine a stress intensity factor. 
For example, the stress intensity coefficient at x ■ 0, when calculated by in- 
tegration, is 1.1413 (table II). This compares with the value 1.122 from which 
the conic section coefficient was calculated. The conic section model is ob- 
viously not self-similar. 

In figure 3(b), the conic-section model and Petroski and Achenbach's model 
provide at best only a rough approximation of the crack profile. Both models 
are essentially two-parameter representations. It seems doubtful that any two- 
parameter representation could closely match the Bueckner or Wigglesworth crack 
profiles. Although more elaborate models could be constructed, the effort may 
not be warranted. Either method may be adequate for most engineering applica- 
tions. The choice of method will nwst likely depend on the information avail- 
able for the reference problem. 


SUMMARY OF RESULTS 

A closed-form expression is given for tho weight function for a strip with 
a single edge crack. The expression is valid for relative crack lenghts from 
zero to unity. Computed weight functions are within 4 percent of reference 
values for relative crack lengths of one-half and less. Computed ‘■tress in- 
tensity coefficients are within 1.8 percent (tension) and 2.7 percent (bending) 
of reference values for relative crack lengths from zero to unity. Crack mouth 
opening displacement coefficients are within 2.1 percent of reference values 
for relative crack lengths of 0.8 and less. 

This paper and (ref. 8) each present methods for developing approximate 
weight functions. Each assumes a form for the crack opening displacements. 
Neither method appears to have an overwhelming advantage, and the choice may 
depend on the problem to be solved and the available information. 


*The special form of the displacement coefficient used in figure 3(b) 
deserves comment. It can be shown that 

(x/a)"^^^ (E'n/oa)|^^g^_Q -^5 Kj/ayira 


Thus the form used in figure 3(b) is particularly useful since the inter- 
cept at x/a « 0 is proportional to the stress intensity factor and the 
intercept at x/a - 1 is the crack mouth displacement. In this way, two 
significant features of the opening shape may be clearly seen on the same 
plot. 
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TABLE 1. - WEI6HT FUNCTIONS FOR THE EDGE-CRACKED STRIP 
( 4 ) Bu«ckn«r (taken from second table in [2] except as noted). 


a/H 

x/a 


0.1 

noon 

0.3 

IQQII 

I^QII 


QQII 

0.8 

0.9 

1.0* 


N(x/a) 

0 . 00 “* 


1.1407 

1.2079 

1.2797 

1.3590 

1.4468 

i.r 24 

1.6433 

1 

1.7453 

1.6424 

.05 

1.0659 

1.1361 

1.2115 

1.2928 


1.4751 

1.5756 

1.6806 

1.7874 

1.9161 

.10 

1.0818 

1.1688 

1.2614 

1.3601 


1.5768 

1.6942 

1.8161 

1.9398 

2.0697 

.15 

1.1061 

1.2176 

1.3348 

1.4581 


1.7248 

1 .3679 

2.0158 

2.1656 

2.3264 

20 

1.139( 

1.2841 

1.4358 

1.5945 

1.7606 

1.9343 

2.1151 

2.3020 

2.4934 

2.6869 

.25 

1.177: 

1.3623 

1.5554 

1.7569 

1.9670 

2.1859 

2.4132 

2.6477 

2.8874 

3.1637 

.30 

1.225! 

1.4610 

1.7077 

1.9653 

2.2338 

2.5132 

2.8033 

3.1032 

3.4107 

3.7650 

35 

1.2831 

1.5814 

1 .8949 

2.2233 

2 . 5665 

2.9424 

3.2962 

3.6811 

4.0772 

4.5178 

.40 

1.3501 

1.7224 

2 1162 

2.5306 

2.9650 

3.4190 

3.8918 

4.3820 

4.8867 

5.4635 

. 4 ; 

1.4348 

1.9024 

2.4012 

2.9296 

3.4864 

4.0705 


5.3160 

5.9737 

6.6653 

.50 

1.5376 1 

2.1237 

2.7550 

3.4288 

4.1430 

4.8956 



7.3649 

8.2149 


Jrhis column calculated from eqn. (2.5) of [1] except for x«0. 
“This row caiculaieti from eqns. (4.32) and (4.35) of [2]. 


(b) This report, Eqs. (3) to (5). 


a/u 

x/a 

0.1 

0.2 

0.3 

0.4 

Di 

mm 


0.8 

0.9 

1.0 

N(x/a) 

0.00 

.05 

.10 

.15 

.20 

.25 

.30 

.35 

.40 

.45 

.50 

1.0768 

1.0836 

1.1004 

1.1242 

1.1452 

1.1911 

1.2352 

1.2882 

1.3521 

1.4300 

1.5263 

1.1548 

1.1686 

1.2022 

1.2505 

1.2924 

1.3875 

1.4793 

1.5906 

1.7262 

1.8933 

2.1019 

1.2343 

1.2549 

1.3056 

1.3787 

1.4417 

1.5894 

1.7322 

1.9069 

2.1218 

2.3885 

2.7242 

1.3152, 

1.3426' 

1.4104 

1.5090 

1.5932 

1.7967 

1.9440 

2.2372 

2.5382 

2.9144 

3.3910 

1.3975 

1.4318 

1.5169 

1.6415 

1.7470 

2.0095 

2.2646 

2.5811 

2.9752 

3.4701 

4.1001 

1.4815 

1.5225 

1.6250 

1.7762 

1.9033 

Z.?V9 

2.5440 

2.9385 

3.4320 

4.0545 

4.8497 

1.5670 

1.6147 

1.7349 

1.9132 

2.0620 

2.4519 

2.8323 

3.3093 

3.9085 

4.6667 

5.6382 

1.6542 

1.7087 

1.8465 

2.0526 

2.2233 

2.6815 

3.1293 

3.6932 

4.4040 

5.3058 

6.4641 

1.7431 

1.8043 

1.9600 

2.1943 

2.3873 

2.9169 

3.4352 

4.0903 

4.9183 

5.9712 

7.3261 



1 .8338 
1.9018 
2.0755 
2.3387 
2.6949 
3.1581 
3.7499 
4.5002 
5.4509 
6.6620 
8.2230 


t^Largest difference, +3.85 pe+cent. 
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TABLE II. - STRESS INTENSITY FACTOR COEFFICIENTS 


a/W 

Unifonn tension, 
(Kj/oK.a) (1 - L)^'^/(l + 2») 

Pure Pending, 
(Kj/o/,a) (1 - 


Ref. [7] 

This report 

Ref. [7] 

This report 

iBfBI 

1.122 

1.1413 

1.172 

1.1413 

.05 

.961 

.9780 

.590 

1.0107 

.10 

.849 

.8631 

.830 

.9122 

,15 

,767 

.7768 

.811 

.8328 

.20 


.7096 

.747 

.7667 

.25 

.653 

.6558 

.6RJ 

.7108 

.TO 

.611 

.6121 

.649 

.6633 

.35 

.577 

.5761 

.610 

.6224 

.40 

.548 

.5457 

.577 

.5869 

,45 

.523 

.5205 

.549 

.5562 

.50 

.501 

.4984 

.523 

.5290 

.55 

.482 

.4790 

.501 

.5047 

.60 

.464 

.4618 

.431 

.4830 

.65 

.449 

.4466 

.463 

.4635 

.70 

.435 

.4329 

.446 

.4461 

.75 

.423 

.4207 

.432 

.4306 

.80 

.411 

.4097 

.418 

.4167 

.85 

.401 

.3999 


.4043 

.90 

.391 

.3908 

.394 

.3932 

.95 

.382 

.3823 

.384 

.3831 


.374 

1 






TABLE III. - CRACK MOUTH 
DISPLACEMENT COEFFICIENTS 
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NORMALIZED DISTANCE FROM CRACK TIP, x/a 


Figure J. - Conic sections. 



RELATIVE CRACK LENGTH, X • aW 
Figure 2 . - Conic section coefficient for tension load. 


ORiG^r 

OF POOR QUALITY 



(b) Special form. 


Figure 3. - Crack opening displacement for 
the edge crack In a half-plane. 


